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Abstract 

The resummation of logarithms in Quantum Field Theories is a 
long tale plenty of successes, yet the resummation of logarithms in 
non-relativistic theories has remained elusive. This was the most frus- 
trating, since the first quantum field theory log ever computed was the 
Lamb shift one. We briefly review recent progress on the resumma- 
tion of logarithms of a, which appear in the physics of non-relativistic 
states, using effective field theories. We put special emphasis on the 
basic formalism. 
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Heavy quark-antiquark systems near threshold are characterized by the 
small relative velocity v of the heavy quarks in their center of mass frame. 
This small parameter produces a hierarchy of widely separated scales: m 
(hard), mv (soft), mv"^ (ultrasoft), ... The factorization between them is effi- 
ciently achieved by using effective field theories, where one can organize the 
calculation as various perturbative expansions on the ratio of the different 
scales, effectively producing an expansion in v. The terms in these series get 
multiplied by parametrically large logs: Inw , which can also be understood as 
the ratio of the different scales appearing in the physical system. Again, effec- 
tive field theories are very efficient in the resummation of these large logs once 
a renormalization group (RG) analysis of them has been performed. There- 
fore, the use of effective field theories has opened the possibility to solve a 
problem, otherwise open, since the first non-relativistic logs appeared around 
fifty years ago (in the Lamb shift in the Hydrogen atom). We will review in 
this paper the recent progress achieved in potential NRQCD (pNRQCD) [T] 
on the issue of the resummation of Inf terms in the weak coupling regime 
(we will obviate any non-perturbative effects in what follows). In particular, 
we will focus on the theoretical aspects. 

Potential NRQCD (pNRQCD) is defined by its particle content and cut- 
off PpNR = {t'p, z^jis}, where Pp is the cut-off of the relative three-momentum of 
the heavy quarks and v^s is the cut-off of the three-momentum of the gluons 
and light quarks. They satisfy the following inequalities: |p| -C z/p ^ m 
and p^/m <^ v^s |p|, where typically |p| ~ mv. Note that no gluons 
or light quarks with momentum of 0(|p|) are kept dynamical in pNRQCD. 
The motivation to integrate out these degrees of freedom is that they do not 
appear as physical (on-shell) states near threshold. Nevertheless, they can 
appear off-shell and, since their momentum is of the order of the relative 
three-momentum of the heavy quarks, integrating them out produces non- 
local terms (potentials) in three-momentum space. Indeed, these potentials 
encode the non-analytical behavior in the transfer momentum of the heavy 
quark, k = p — p', of the order of the relative three-momentum of the heavy 
quarks. 

In this paper, we will mainly consider the situation v ^ as- It should be 
clear however that pNRQCD (in the weak coupling regime) is also valid in 
the situation mv'^ ^ ma^ as long as v is an small parameter. Nevertheless, 
the Coulomb resummation is not necessary in this case. 

Formally, we can write the pNRQCD Lagrangian as an expansion in 



1 



l/r(= 1/r, p) and 1/m in the following way: 



n=-l "^n^2 " 

where the above operators and matching coefficients should be understood 
as bare. As for the renormalized quantities, we define V as the potentials 
and V as the (almost) dimensionless constants in it. The latter are in charge 
of absorbing the divergences of the effective field theory. Therefore, they will 
depend on Up and z/^^. At next-to- leading order in the multipole expansion 
the Lagrangian can be written as^ 

LpNRQCD = J d\d'Rtri^S^{tdo - K}s + 0^{tDo - K}0 (2) 
+gVAr ■ OE(R, t)S^ + qVat ■ O^E{R, t)S 
+ 1 V^r • 00tE(R, t) + ^Vbt ■ O^OE{R, t) | , 

where hg and ho are the singlet and octet quantum mechanical Hamiltonian. 
For illustration, at low orders (see Ref. Pj for notation and further details), 
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in the equal mass case, where Cf = {N^ — l)/{2Nc) and we will set = C4 = 
1. A similar expression holds for the octet Hamiltonian changing the label 
of the matching coefficients {s o). 

The matching process, which basically means the computation of the 
potentials, is carried out for a given external incoming (outcoming) momen- 
tum p (p'). Therefore, one has to sum over all of them in the pNRQCD 
Lagrangian, since they are still physical degrees of freedom as far as their 



^iS* and O stand for the heavy quarkonium bihnear field in a singlet and octet config- 
uration under ultrasoft gauge transformations. They both depend on the relative, r, and 
center of mass, R, coordinate. 
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momentum is below Up. In position space, this means that an integral over 
X, the relative distance between the heavy quarks, appears in the Lagrangian 
when written in terms of the heavy quark-antiquark bilinear field. 

Within pNRQCD, integrals over p (or x) appear when solving the Schro- 
dinger equation that dictates the dynamics of the heavy quarkonium near 
threshold. At low orders, these integrals are finite effectively replacing p by 
~ mas- Nevertheless, at higher orders in quantum mechanics perturbation 
theory and/or if some singular enough operators are introduced (as it is the 
case of the heavy quarkonium production currents) singularities proportional 
to In Vp appear. These must be absorbed by the potentials or by the matching 
coefficients of the currents. 

Let us now describe the matching between QCD and pNRQCD within 
an RG framework. We first address the procedure that gives the running 
of the potentials. One first does the matching from QCD to NRQCD. The 
latter depends on some matching coefficients: c(z/s) and divp, Us), which can 
be obtained order by order in a<j (with z/p = Us) following the procedure 
described in Ref. for the c(z/s) and j3j for the d^Up, Vg)- is the ultraviolet 
cutoff of the three-momentum of the gluons in NRQCD. The c{y^ stand for 
the coefficients of the operators that already exist in the theory with only one 
heavy quark (ie. HQET) and the (i(z/p, stand for the coefficients of the four 
heavy fermion operators. The starting point of the renormalization group 
equation can be obtained from these calculations by setting Vp = = m 
(up to a constant of order one). In principle, we should now compute the 
running of Vp and v^. The running of the civ^ can be obtained using HQET 
techniques [S]. The running of the divp^ Ug) is more complicated. At one-loop, 
Up does not appear and we effectively have di^Up, Vg) — d{vs), whose running 
can also be obtained using HQET-like techniques [2.. At higher orders, the 
dependence on Up appears and the running of the d^Up, Vg) becomes more 
complicated. Fortunately, we need not compute the running of d in this 
more general case because, as we will see, the relevant running of the d for 
near threshold observables can be obtained within pNRQCD. 

The next step is the matching from NRQCD to pNRQCD. The latter 
depends on some matching coefficients (potentials). They typically have 
the following structure: V^(c(i/s), (i(i/p, i^s), z/^, z/„s, r). After matching, any 
dependence on i/g disappears since the potentials have to be independent of Vg- 
Therefore, they could be formally written as V^(c(l/r), d{vp, 1/r), 1/r, z/^^, r). 
These potentials can be obtained order by order in Ug following the procedure 
of Refs. PI E] . The integrals in the matching calculation would depend on 
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a factorization scale /z, which should correspond either to or to u^s- In 
the explicit calculation, they could be distinguished by knowing the UV and 
infrared (IR) behavior of the diagrams: UV divergences are proportional to 
In i/g, which should be such as to cancel the Ug scale dependence inherited from 
the NRQCD matching coefficients, and IR divergences to Uus- In practice, 
however, as far as we only want to perform a matching calculation at some 
given scale fJ^ = i^s = ^us, it is not necessary to distinguish between and Uus 
(or if working order by order in ag without attempting any log resummation) . 

Before going into the rigorous procedure to obtain the RG equations of 
the potentials, let us first discuss their structure on physical grounds. The 
potential is independent of Ug. This allows us to fix z/^ to 1/r that, in a way, 
could be understood as the matching scale for Ug"^. Therefore, 1/r, the point 
where the multipole expansion starts, would also provide with the starting 
point of the renormalization group evolution of Uus (up to a constant of 
order one). The running of can then be obtained following the procedure 
described in Refs. j71l2]. Formally, the renormalization group equations of 
the matching coefficients due to the z/^j^-dependence read 

iyus-^V = By{V). (4) 

From a practical point of view one can organize the RG equations within 
an expansion in 1/m. At 0(l/m°), the analysis corresponds to the study of 
the static limit of pNRQCD, which has been carried out in Ref. 0. Since 
V-i 7^ 0, there are relevant operators (super-renormalizable terms) in the La- 
grangian and the US RG equations lose the triangular structure that we en- 
joyed for the RG equations of Ug. Still, if V-i <^ 1, a perturbative calculation 
of the renormalization group equations can be achieved as a double expansion 
in V-i and Vq, where the latter corresponds to the marginal operators (renor- 
malizable interactions)). At short distances (1/r ^ Aqcd), the static limit 
of pNRQCD lives in this situation. Specifically, we have V^i = {av^,avo}j 
that fulfills V-i ~ Q:g{r) -C 1; Vq = c^si^us) and Vi = {Va,Vb} ~ 1. There- 
fore, we can calculate the anomalous dimensions order by order in agluus). 
In addition, we also have an expansion in Vli. Moreover, the specific form 
of the pNRQCD Lagrangian severely constrains the RG equations general 

^In practice, the potential is often first obtained in momentum space so that one could 
then set = k. Note, however, that this is not equivalent to fix = 1/r, since finite 
pieces will appear after performing the Fourier transform. 
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x> 



i/(E - - pVm) 

Figure 1: The UV divergences of this diagram in pNRQCD give the leading 
non-trivial ultrasoft running of a^^, -D^^-*, D'f'\ 



structure. Therefore, for instance, the leading non-trivial RG equation of 

tty^ reads 

^^s^^^v. = l^Vl ((^ - ay^ + C,ay}j' + OiVt.V,, V^Vl,) . (5) 

At higher orders in 1/m the same considerations than for the static limit 
apply here as far as the non-triangularity of the RG equations is concerned. 
In general, we have the structure (V^^-* = Vm) 

^us-^y^^- E with ^n, = n,E^.=^, 

'^^"s {ni}{'mi} i=l i=l 

(6) 

and one has to pick up the leading contributions from all the possible terms. 
Actually, as far as the NNLL heavy quarkonium mass is concerned, the rele- 
vant ultrasoft running can be obtained by computing the diagram displayed 
in Fig. n (see |BJ (one also has to consider the running of Va, which happens 
to be zero). 

At the end of the day, we would have V^(c(l/r), ^(t'p, 1/r), 1/r, z/„s, r), 
where the running on v^s is known and also the running in 1/r if the d is 
z/p- independent. So far, the only explicit dependence of the potential on Vp 
appears in the d. Nevertheless, the potential is also implicitly dependent on 
the three-momentum of the heavy quarks through the requirement 1/r ~ 
p -C z/p, and also through z/^,s, since Vus needs to fulfill p^/m <C i/^s <^ |p| 
in order to ensure that only soft degrees of freedom have been integrated 
out for a given |p|. This latter requirement holds if we fix v^s = Vp/m (this 
constraint tells you how much you can run down v^s in the potential before 
finding the cutoff z/^ /m caused by the cutoff of p) . 
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Within pNRQCD, the potentials should be introduced in the Schrodinger 
equation. This means that integrals over the relative three-momentum of 
the heavy quarks take place. When these integrals are finite one has p ~ 
1/r ~ mas and p^/m ~ maf. Therefore, one can lower Uus down to ~ ma^ 
reproducing the results obtained in Ref. [2]. In some cases, in particular in 
heavy quarkonium creation, the integrals over p are divergent, and the log 
structure is dictated by the ultraviolet behavior of p and 1/r. This means 
that we can not replace 1/r and by their physical expectation values but 
rather by their cutoffs within the integral over p. Therefore, for the RG 
equation of Up, the anomalous dimensions will depend (at leading order) on 
V^ci^Up), diup, Up), Up, Up/m, UpY and the running will go from Up m down 
to Up ~ mas- Note that, at this stage, a single cutoff, Up, exists and the 
correlation of cutoffs can be seen. The importance of the idea that the 
cutoffs of the non-relativistic effective theory should be correlated was first 
realized in Ref. p. 

With the above discussion in mind, the matching between NRQCD and 
pNRQCD could be thought as follows. One does the matching by computing 
the potentials order by order in at the matching scale Up = u^ = u^s 
following the procedure of Refs. [H IH] (by doing the matching at a generic 
Up some of the running is trivially obtained). The structure of the potential 
at this stage then reads V{c{up),d{up,Up),Up,Up,Up) (and similarly for the 
derivatives with respect Inr of the potential). This provides the starting point 
of the renormalization group evolution of u^s (up to a constant of order one). 
The running of u^s can then be obtained following the procedure described in 
Refs. [71I2]. For the final point of the evolution of Uus, we choose Uus = Vp/m. 
At the end of the day, we obtain V{c{up), d{up, Up), Up, u'^/m, Up) = V{up). 

The running of Up goes from Up = m (this was fixed when the matching 

■^Roughly speaking, this result can be thought as expanding Inr around Inup in the 
potential ie. 



V{c{l/r),d{iyp,l/r),l/r,iyp/m,r) ~ V{c{vp),d{vp,Vp),Up,vl/m,Vp) 

+ •■•■ (7) 



+ \n{vpr)r—V 
ar 



The ln(j/pr) terms give subleading contributions to the anomalous dimension when intro- 
duced in divergent integrals over p. A more precise discussion would require a full detailed 
study within dimensional regularization. An explicit example of this type of corrections 
appear in the computation of the hyperfine splitting of the heavy quarkonium at NLL 

[iniini. 



6 



between QCD and NRQCD was done) up to the physical scale of the problem 
Up ~ mag. If the running of the NRQCD matching coefficients is known, the 
above result gives the complete running of the potentials. The procedure 
to get the running of the c is known at any finite order. For the d it is 
just known at one-loop order, since, at this order, it is only i/g-dependent. 
Nevertheless, at higher orders, dependence on i/p appears. Therefore, the 
above method is not complete unless an equation for the running of Up is 
provided. This is naturally given within pNRQCD. It appears through the 
iteration of potentials. Let us consider this situation more in detail. The 
propagator of the singlet is (formally) 

(8) 



At leading order (within an strict expansion in ag) the propagator of the 
singlet reads 

= Gc{E) 



E - E - p7m - Cfas/r ' 

If we were interested in computing the spectrum at 0(mQ;f), one should 
consider the iteration of subleading potentials (Shg) in the propagator as 
follows: 

G,iE)ShgG,{E)---5hsG,iE). (9) 

In general, if these potentials are singular enough, these contributions will 
produce logarithmic divergences due to potential loops. These divergences 
can be absorbed in the matching coefficients, D^^] and -D^?^, of the local 
potentials (those proportional to the (5(^^(r)) providing with the renormal- 
ization group equations of these matching coefficients in terms of i^p. Let us 
explain how it works in detail. Since the singular behavior of the potential 
loops appears for |p| ^ «s/r, a perturbative expansion in as is allowed in 
Gc{E), which can be approximated by 

I I ^G(^\E) ^ 



E - p'^/m 



Therefore, a practical simplification follows from the fact that the Coulomb 
potential, —Cf^, can be considered to be a perturbation as far as the com- 
putation of the In i/p ultraviolet divergences is concerned. This means that 
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the computation of the anomalous dimension can be organized within an 
expansion in and using the free propagators G^^\ Moreover, each G^^^ 
produces a potential loop and one extra power of m in the numerator, which 
kills the powers in 1/m of the different potentials. This allows the mixing of 
potentials with different powers in 1/m. One typical example would be the 
diagram in Fig. 2. The computation of this diagram would go as follows: 



(E) ^J/'^ 6^'^ (r)GW {E)Cf^G^^^ (E) 6^'^ (r)^^ (E) . (10) 



Using 5(3) (r) = |r = 0)(r = 0|, we can see that the relevant computation 
reads (instead of ay^ one could use since the non-trivial running of ay^ is 
a subleading effect. Nevertheless, we keep ay, since it allows to keep track 
of the contributions due to the Coulomb potentials) 

(r = 0|Gf(E)C;^Gf(i?)|r = 0) (11) 
d'^p' f d'^p m inav m m^av 1 



2tiY J {27tY p'^ - niE q2 p2-mE ^ Wn e' 

(2) 

where D = 4 + 2e and q = p — p'. This divergence is absorbed in ^ 
contributing to its running at next-to-leading-log (NLL) order as follows 

^p^^D?i(^p) ~ »vMD?i\i^p) + ■■■■ (12) 

Therefore, even without knowing the running of the d (which need to be 
known at NLL order in this case), we can obtain the running of the potential 
(one can also think of trading Eq. (fT^ into an equation for d, which is the 

(2) 

only unknown parameter within the potential). This is so because -D)^^ is 
only needed with LL accuracy in the right-hand side of Eq. p2|). 

This finishes the procedure to the RG equations. The above method deals 
with the resummation of logs due to the hard, soft and ultrasoft scales. Nev- 
ertheless, for some specific kinematical situations even smaller scales could 
appear. Their study, however, has not yet been carried out. In any case, 
pNRQCD can be considered to be the right starting point to study these 
kinematical situations. 

This line of investigation has lead to several new results on heavy quarko- 
nium physics. They can be summarized in the following way 
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d.s Vs d,s 



Figure 2: One possible contribution to the running of D)^ l at NLL. The 
first picture represents the calculation in terms of the free quark- antiquark 
propagator G^^^ and the potentials (the small rectangles) . The picture below 
is the representation within a more standard diagrammatic interpretation in 
terms of quarks and antiquarks. The delta potentials are displayed as local 
interactions and the Coulomb potential as an extended in space (but not in 
time) object. 



• The correction to the heavy quarkonium energy at NNLL [2], i.e. cor- 
rections of order 

6E ~ ma^ + ma^ In a + ma^ In^ a + ■ ■ • (13) 

• Corrections to the heavy quarkonium hyperfine sphtting at LL |2] (first 
obtained in Ref. |I2]) and NLL HUIIII] 

5Ehf ~ ma^ + ma^ \n a + ma^ \yi' a + ■ ■ ■ (14) 
+ ma + ma In a + ma In a + ma In a + ■ ■ ■ 

• The decays are known with NLL accuracy (this result can be easily 
applied to t — t production threshold or non-relativistic sum rules since 
the running of the electromagnetic current matching coefficient is the 
only non-trivial object that appears at NLL running) |13j 

T{VQ{nS) ^ e-^e-) ~ ma^(l + In a + In^ a + ■ ■ ■) (15) 
r(PQ(n5') 77) ~ ma^(l + a^lna + a^ln^a H ) 
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and for the ratio with NNLL accuracy ^3] 

~ 1 + a^lna + a^ln^a H (16) 



TiyginS) ^ e+e- J ^ , , 2 



T{PQ{nS) 77) 

In a + In^ a + 



There has also been a lot of work on the resummation of logarithms using 
the vNRQCD framework. Unfortunately, its first formulation suffered from 
some mistakes (in particular concerning the treatment of ultrasoft modes), 
which lead to incorrect results for the heavy quarkonium mass at NNLL ^2] 
and the electromagnetic current matching coefficient at NLL Fortu- 
nately, they have been solved in Ref. jTH] and their results now agree with 
those obtained in pNRQCD [21 US]- 

These results may have an important phenomenological impact in sev- 
eral situations. Let us enumerate a few of them. The determination of the 
bottom and charm masses (using the experimental value of the ground state 
heavy quarkonium masses or non-relativistic sum rules). The determination 
of the rib{lS) mass, the hyperfine splitting of the ground state Be system, 
or theoretical improved determinations of the rjc- One can also try to ob- 
tain improved determinations for the inclusive electromagnetic decays of the 
heavy quarkonium. On the other hand the application of this program to t-i 
production near threshold at the Next Linear Collider is one of the main 
motivations to undergo these computations. Moreover, it would also be in- 
teresting to try to perform the resummation of logarithms in semi-inclusive 
radiative decays [T7j. A review on the phenomenological consequences of 
these results is presented elsewhere. 
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